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Photometric variability in FU Ori and Z CMa as observed 
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ABSTRACT 

Photometric observations obtained by the MOST satellite were used to characterize 
optical small scale variability of the young stars FU Ori and Z CMa. Wavelet analy- 
sis for FU Ori reveals the possible existence of several 2-9 d quasi-periodic features 
occurring nearly simultaneously; they may be interpreted as plasma parcels or other 
localized disc heterogeneities revolving at different Keplerian radii in the accretion 
disc. Their periods may shorten slowly which may be due to spiralling in of individual 
parcels toward the inner disc radius, estimated at 4.8 ± 0.2 R Q . Analysis of additional 
multicolour data confirms the previously obtained relation between variations in the 
B — V colour index and the V magnitude. In contrast to the FU Ori results, the oscil- 
lation spectrum of Z CMa does not reveal any periodicities with the wavelet spectrum 
possibly dominated by outburst of the Herbig Be component. 

Key words: star: individual: FU Ori, Z CMa, stars: accretion: accretion discs. 



1 INTRODUCTION 

This paper presents results obtained for two well known 
young stars, FU Ori and Z CMa, as a continuation of the 
MOST satellite photometric variability studies of Young 
Stellar Objects. 

FU Ori is an object of special interest which has been 
known since 1937, when in a timescale of one year its bright- 
ness rose from 16 to 9.5 mag (in the photographic blue-band 
system) and t hen started to d e cay slowly at the r ate of 
0.015 mag/yr <|Wachmannlll939l ; iKenvon et alJl200Ch . This 
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outburst, and two similar events observed in the early 1970's 
in V1057 Cyg and V1515 Cyg, led to t he creation of a class 
of eruptive young stars - "FUor" stars dHerbid 19771) - cur- 
rently cons i sting of about 20 members 1 Semkov fc Peneval 
12011 12012L iMiller et aill201l| j. The pre-outurst spectra of 
FUors obtained for two members of the class, V1057 Cyg 
and V2493 Cyg (HBC 722), revealed that the progenitors 
were Classical T Tauri- type stars (CTTS). 

An enhanced accretion in the disc at a rate of 
about 10~ 4 M@/yr, resulting in a major increase of the 
disc surface brightness and dominating over the stellar 
flux, was proposed as the s ource of the FUor outbursts 
l|Hartmann fc Kenvonl Il985l . 1 19961 ). Disc overy of compan- 
ions to FU Ori, Z CMa and other FUors (|Wang et al.ll2004l : 
iReipurth fc Aspin|[200l ). has initiated discussion whether 
FUor outbursts may be triggered by pert urbations in the ac- 
cretion discs at close periastro n passage (|Bonnell fc Bastienl 
Il992l : IReipurth fc Aspinll2004l ) . 

During an outburst of a FUor, emission lines, which are typi- 
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cal for CTTSs, almost completely disappear while the visual 
spectrum is dominated by absorption features produced in 
an inner accretion disc radiating as a stellar atmosphere of 
a F-G supergiant star. The outer, colder parts of the FUor 
discs produce a K-M type supergi ant spectrum observable 
in the infrared (jKenvon et al.l !l988l In accordance with the 
location of the brightest parts of the Keplerian disc seen 
in different wavelengths, the absorption line broadening di- 
minishes from the visual to the near- and mid-infrared spec- 
tral re g ions, as observed f or FU Ori by_ Kenvon fc Hartmannl 
|l989MHartmann et all (|2004h and lZhu et all i|2009l l 



Z CMa was recognized as a young star by Herbig 

(1960). The discovery, using infrared speckle interferome- 
try, that the star is a very close 0.1 arcsec visual binary 
<|Koresko et al.lll99ll ; lHass et al.lll993l l resolved the initial 
difficulties in inter pretation of the co mplex spectral prop- 
erties of the star (|Covino et alj|l984f ). Spectropolarimeric 
observations revealed that it is the Herbig Be component 
of the binary which dominates the infrared continuum and 
total luminosity of the system. The same star is appar- 
ently the source of the emission lines polarized in the dusty 
disc envelope and ob served in visual part of the spectrum 
I Whitney et all 1 1993( 1 The Herbig B e star is also responsi- 
ble for AV = 1 — 3 mag outburs ts (|van den Ancker et all 
|2004| ; iGrankin fc Artemenkoll2009ll caused by variable scat- 
tering geometry (ISzeifert et alll2010f). pe riods of strong mass 
loss fr om the disc (|Benistv et all 120101 ). or the EXor out- 
burst l|Whelan et all l2010l l These outbursts are superim- 
posed on the slowly decaying light curve of the visually 
brighter FUor star. This component contributes a typical 
- for FUors - rotationally broa dened absorption s pectrum 
produced in the accretion disc (|Weltv et alj 1 19921 1 it sup- 
plied about 80% of the observed unpolarized flux at visual 
wavelengths (|Whitnev et alj|l993h when the system bright- 
ness w^s_K_=_9 ; 7mag_ ; _^ 

iKenvon et all |2000f) argued that similarly to cata- 
clysmic variables, accretion discs of FUors should produce 
nickering variability as a characteristic signature observable 
at visual wavelengths. By comparing a limited amount of 
ground-based data to the Monte-Carlo synthetic variability 
model the authors suggested the possibility of ~ 1 d quasi- 
periodicities with a U-band amplitude of about 0.035 mag, 
which would originate in the inner edge of the FU Ori 
disc. Motivated by these results and by the lack of similar 
studies for Z CMa, we decided to re-examine this issue by 
means of continuous, high precision, space-based photomet- 
ric observations. Additionally, to obtain information on the 
wavelength- amplitude dependence of the flickering, we ob- 
served FU Ori simultaneously by means of a ground-based 
telescope using intermediate-width Stromgren v and b fil- 
ters. We describe details of these observations in Section [21 
The methods used for light curve analysis and the results 
obtained for our targets are presented in Section [3] and then 
summarized in Section [5] 



2 OBSERVATIONS AND DATA REDUCTIONS 

The optical system of the MOST satellite consists of a 
Rumak-Maksutov f/6, 15 cm reflecting telescope. The cus- 
tom broad-band filter covers the spectral range of 380 - 
700 nm with the effective wavelength falling close to the 



Johnson V band. The pre-launch ch a racter istics of the mis- 
sion are described by I Walker et all f2 0031 and the initial 
post-launch performance bv lMatthews et all (|2004l l 

The stars investigated in this paper were observed in the 
direct-imaging mode of the satellite. FU Ori was observed 
nearly continuously for 28 days between 13 December, 2010 
and 9th January, 2011, during 362 satellite orbits. The indi- 
vidual exposures were 30 s during the first part of the run 
and 60 s during the last 10 days of the run. Immediately 
after the FU Ori observations, on 10 January 2011 MOST 
started a 13 day long monitoring of Z CMa, which was ob- 
served with 60 s exposures. Some occasional interruptions 
in data acquisition, visible in the light curves (Figure [T]) did 
not impact the scientific results. 

The dark and flat calibration frames for the MOST 
data were obtained by averaging a dozen empty-field images 
specifically taken during each observing run, or - for the case 
of the 30 sec long exposures of FU Ori - from frames with the 
target localized far beyond its optimal position due to occa- 
sional satellite guiding errors. Aperture photometry of the 
stars was obtained from the dark and flat corrected images 
by means of the DAOPHOT II package (|Stetsonlll987T l. As in 
our previous investigations, a weak correlation between the 
star flux and the sky background level within each MOST 
orbit was noted and removed; it was most probably caused 
by a small photometri c nonlinearity in the electronic system 
fsee lSiwak et al.ll2010l ). As a result, we obtained very good 
quality light curves, particularly for Z CMa (see Fig. [T]). 
In the analysis we used the satellite-orbit (101 min) aver- 
ages for FU Ori whose median error was 0.0036 of the mean 
normalized flux. We note that FU Ori is surrounded by a 
bright nebula which slightly decreased the photometric ac- 
curacy of the data. For the brighter Z CMa, the data were 
binned into smaller, 14.5 min bins consisting of typically 10 
- 15 observations with the median error 0.0011 per bin. The 
finer binning was done in order to preserve information on 
the short time-scale variability of the star. 

During four nights of Jan 4, 5, 8, 9, 2011, we simultane- 
ously observed FU Ori in v and b Stromgren filters using the 
60 cm telescope at the Mount Suhora Observatory, Cracow 
Pedagogical University. We were unable to obtain similar 
data for Z CMa due to poor weather conditions. The data 
were reduced in a standard way in the ESO-MIDAS soft- 
ware environment. Aperture photometry was obtained with 
the DAOPHOT II package. The photometry of FU Ori was 
done differentially utilizing a mean comparison star made 
of the three nearby, somewhat fainter stars having colour 
indices similar to FU Ori (Tab. [T]). Thanks to the proxim- 
ity of the stars and similarity of the colours, corrections for 
differential and colour extinction were unnecessary. 



3 RESULTS OF THE LIGHT-CURVE 
ANALYSIS 

W e performed an a lysis of the MOS T data in a s imilar way 
to iRucinski et all (20081 . l201Gv i and ISiwak etall (|201ll l i.e. 
the Fourier analysis was done by simple least-squares fits 
of expressions of the form /(/) = c o(f) + c\(f)cos[2ir(t — 
to)f] + C2(f) sin[2ir(t — to) f] for an appropriate range of fre- 
quencies /, with the step of Af — 0.01. The amplitude a(f) 
for each frequency was evaluated as the modulus of the pe- 
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Figure 2. Fourier analysis of the variability data computed from the mean-orbital data points of FU Ori (left panel) and from all data 
points of Z CMa (right panel). The amplitude errors estimated through bootstrap repeated sampling are represented by small points. 
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Figure 3. Morlet-6 wavelet spectra of FU Ori (left panel) and Z CMa (right panel) calculated for the whole accessible period range up 
to 11 and 5 days, respectively. The light curve in normalized flux units is shown at the bottoms of the panels. Ranges not affected by 
edge effects in the wavelt transformation are located between the two white broken lines. 
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Table 1. Basic informations on FU Ori and the comparision stars 
used for differential photometry. 

*The average V magnitude of FU Ori during MOST observa- 
tions as obtained from observations made by Konstatin Grankin 
(priv.comm.), Bruno Alain and Timar Andras (AAVSO). 



star 


V 


B-V 


FU Ori 


9.74* 


1.279(47) 


GSC 00714-0203 


10.555(88) 


1.069(85) 


GSC 00715-0188 


10.682(93) 


1.646(241) 


GSC 00715-0123 


11.337(119) 


0.960(194) 



riodic component, a(f) = yc[(f) + cfCf)- The bootstrap 
sampling technique permitted evaluation of mean standard 
errors of the amplitudes from the spread of the coefficients 
a,i. This technique, for a uniform temporal sampling - as in 
our case - may give too pessimistic estimates of errors, but 
consistently we prefer this conservative approach. 

For the Fourier analysis of the FU Ori variability, we 
used 362 mean-orbital data points, while in the case of 
Z CMa, which is much brighter, all 15404 measurements, 
obtained every 60 sec were used. For the wavelet analysis 
of Z CMa, the data have been averaged and partly interpo- 
lated into 1291 points on an equal-spacing grid at 0.010067 d 
(14.5 min) steps. Similar interpolation for FU Ori, which 
suffered more interruptions, led to mapping into a grid of 
395 points. The Morlet-6 wavelet provided the best match 
between the time-integrated power spectrum and original 
frequency spectrum of both FU Ori and Z CMa. This was 
noticed by IRucinski et all (120081) and encoun tered later by 
iRucinski et all l|201Ch and ISiwak et all (|201ll ). The Morlet 
transforms of other orders result in systematic differences in 
the period scale. 



3.1 FU Ori 

3.1.1 Analysis of the Fourier and the wavelet spectrum 

The upper envelope of peak amplitudes in the Fourier trans- 
form of the FU Ori data (Fig. [2l left panel) appears to scale 
as flicker noise, oc l/^/J, where / is the frequency. The first 
most prominent maximum at / = 0.107 c/d is visible as the 
largest and most intense feature in the wavelet spectrum 
(Fig. [3] left panel) and directly in the light curve. Its period 
shortens gently from ~9 to ~8 days, somewhat similarly 
to th e case of TW Hya (jRucinski et all |200S| ; ISiwak et all 
2011). However, we note that this periodicity is based on 
only three cycles so that its reality may be questioned. To 
test the validity of our results, we conducted tests on re- 
covery of drifting, quasi-periodic features which additionally 
change their amplitudes. In each case we obtained very sat- 
isfactory reproduction of synthesised light curves, including 
the rates of period and amplitude changes. 
Following our interpretation of TW Hya, we assume that 
the oscillation feature is produced by plasma condensations 
and/or other accretion disc heterogeneities, in turn produced 
by interactions of stellar magnetic field with the inner disc 
plasma, although we have no clear picture what causes these 
instabilities. If this is a correct view, then the period vari- 
ations reflect a change in localization within the disc at 
radii (estimated from Keplerian periods) from about 12.2 to 



11.3 Rn, assuming a central star mass of 0.3 Mq (|Zhu et al.l 
2007) . The low incl ination of the accretion disc of 55 deg 
(Mal beTet al.| [2005) should indeed permit to observe such 
inner disc structural changes. 

Other peaks in the Fourier spectrum (Figure [2} local- 
ized at f = 0.166, 0.215 and 0.263 c/d, (P = 6.0, 4.7 and 
3.8 d, respectively) are represented by the much fainter and 
apparently constant features visible in different parts of the 
wavelet spectrum. They cannot be well characterized and 
may actually represent processing artifacts for the relatively 
short dataset or a pure random flickering. We note that these 
three periodicities are not directly visible in the FU Ori light 
curve. 

In addition to the features described above, the light 
curve contains a well defined directly visible short-periodic 
signal at / = 0.424 c/d, which disappears after HJD ~ 
2, 455, 560. To improve its visibility we prepared a second 
wavelet image (FigHJ), which reveals one apparently drifting 
feature with two well-defined peaks: the first one is local- 
ized at P = 2.4 d and the second at P = 2.2 d. Within the 
model of orbiting plasma in the inner disc this would cor- 
respond to the Keplerian radii within 5.1 to 4.8 Kq. We do 
not see any obvious short-period signals at periods shorter 
than P = 2.2 d which may indicate that the inner accretion 
disc is truncated at about 4.8 Rq. This value is in agree- 
ment with the inner disc radius of 5.5j^ '« Rr o btain ed from 
interferometric observations by iMalbet et al.l (120051) but is 
3 time s smaller from that obtained by Eisner fc Hillenbrand! 
l|201ll) - 15 ±4 R . 



3.1.2 Analysis of multicolour observations 

Although the weather conditions during the ground-based 
observations were poor, we were very fortunate to have the 
four clear nights exactly when the star showed both the max- 
imum and the minimum of its 8 d long oscillation (Fig[SJ). 
This enabled us to address the wavelength- amplitude depen- 
dence in the b (4670A) and v (4110A) filters of the Stromgren 

photometry. 

As discussed bv lKenvon et al l (|2000h . the observed vari- 
ability amplitudes decrease at shorter wavelengths. In or- 
der to compare our results w ith the relation of th e changes 
A(B - V) and AV found bv lKenvon et all l|2000l ). we used 
the average magnitudes in our two filters and found the 
mean colour index A(0.5x [v + b])- AV = 0.016±0.002 mag. 
It i s larger by 0.06 mag than that returned by Equation 6 of 
the lKenvon et al.1 (2000) paper, but agrees with their finding 
that the B — V colour index becomes redder as the star be- 
comes brighter. This relation allowed the authors constrain 
the effective spectrum producing the flickering to F7-G3 
spectral type. 

The 0.06 mag difference in the colour index can be ex- 
plained by the difference in the location of the dominant 
source of the observed stellar flux. The respective annuli in 
the disc, at different Keplerian distances, should produce 
different wavelength- amplitude relations. The 9-8 d oscil- 
lation feature is assumed to be produced at the distance 
of 12 — 11 Rp,, where according to FU Ori disc models 
|Zhu et al.ll2007l ). T e ff « 1500 K, so that the maximum 
energy is located at A w W/im. Obviously, the amplitudes 
of flux variations arising at such a large distance are con- 
siderably reduced by the dominant flux produced in the in- 
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HJD - 2,455,500 

Figure 4. The upper panel shows the wavelet spectrum of FU Ori 
(as in Fig. |3j but limited to 3 days to improve visibility of shorter 
oscillations. The light curve in normalized flux units is shown at 
the bottom of the panel. Ranges not affected by edge effects in the 
wavelt transformation are located between the two white broken 
lines. 



nermost disc. As the feature approaches the central star on 
a spiral orbit, it moves to warmer disc annuli. As a conse- 
quence, the wide-band MOST filter receives initially a small 
fraction of the energy emitted at 12 R©, but then this frac- 
tion increases as the spiralling-in process continues. 

Changes in the amplitude of the 9-8 day periodic fea- 
ture may be questioned in view of the duration of the whole 
MOST observations lasting 4 weeks. However, the general 
trend of the progressively growing amplitude as the period 
shortens seem to be visible both, directly in the light curve 
(Fig [1] left panel) and in the wavelet spectrum (Fig [3] left 
panel). The amplitude reached about 0.03 mag in the first 
and the second cycle (the overlapping non-sinusoidal os- 
cillations do not allow for a more accurate estimate), and 
0.075 mag in the third. Once again we stress, that this find- 
ing bases on three oscillations only. However, our data do 
reveal a similar effect, of the amplitude growing up from 
0.005 mag to 0.02 mag, also for the P = 2.4 - 2.2 d 
feature; such a drift would correspond to a motion from 
5.1 to 4.8 Rr, where the disc annuli have temperatures 
T cff w 3400 K dZhu et al.ll2007ri. T he difference in the am- 
plitudes with lKenvon et al.l ( 2000T ) may have resulted from 
averaging of independent wave trains during their observa- 
tions or may reflect a real, physical difference in sizes of the 
gas elements. 



3.2 Z CMa 

As mentioned in Section[TJ about 80% of the unpolarized vi- 
sual flux of Z CMa is produced by the FUor component. The 
direct flux from the Herbig Be star is obscured by a dusty 
envelope and visible as scattered continuum light, which is 
polar ized at about 6% and as emission lines (|Whitnev et al.l 
1993). The MOST satellite does not have any polarization 
capabilities and measures the total flux. The observations in 
January 2011 were obtained during the light maximum of a 
new outburst when, according to the AAVSO database, the 
average total system brightness was V — 8.7 mag. This may 
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Figure 5. A fragment of FU Ori light curve covered by the 
ground-based data. The light curves from MOST (small crosses) 
and in b and v filters (squares and circles) were shifted to the 
same level at the light minimum for an easier comparison. 



explain the chaotic shape of the wavelet transform (Fig. [3] 
right panel) containing many incoherent variations down to 
time scales of about 0.3 day. At least a part of these os- 
cillations may have their origin in the night-to-night and 
hour-to-hour variations of H a , H^, N a I D P Cygni absorp - 
tions and in the H/3 emission peak l|Chochol et al.lll998l ). 
Also, from speckle observatio ns obtained in the infrared 
H&cK filters lHass et"al] l|l993h reported that both compo- 
nents of the close visual double were varying independently 
by ~ 0.5 mag. 

The time variability spectrum of Z CMa is some- 
what similar to tha t of the Herbig Ae star HD 37806 
l|Rucinski et alj|20ich . The chaotic variability at all time 
scales and lack of any regularity in the Fourier or wavelet 
transforms strongly suggest that the time variability spec- 
trum is dominated by the Be component which contributes 
at least 20% to the total flux. Lack of any firm results is com- 
pounded by the relatively short duration of the observation 
run which lasted only 13 days. Unfortunately, the MOST 
observations of Z CMa done one year earlier, in 2010, ob- 
tained in a time of relatively moderate light variations were 
of a very poor quality due to technical problems. 



4 SUMMARY 

Continuous MOST satellite observations of FU Ori and 
Z CMa confirmed presence of intensive, rapid (time scales of 
single days) variations of their light which may have general 
characteristics of the flicker noise with amplitudes scaling 
as: a oc 1/vT, similar to the Classical T Tauri star TW Hya 
ijRucinski et al.ll2008l : ISiwak et all 1201 if l. For FU Ori, the 
Fourier and wavelet transform spectra show quasi-periodic 
features, one well defined and possibly changing its period 
within 9 to 8 day, and one much fainter, at about 2.4 days 
drifting down to 2.2 days. But we note an a bsence of any 
1-day periodicity in our data, contrary to the lKenvon et aD 
(2000) results which were obtained from a limited amount 
of unequally-spaced, poor-quality data. We tentatively in- 
terpret these quasi-periodic variations as produced by hot 
plasma condensations and/or disc heterogeneities which de- 
velop in the magneto-rotationaly unstable inner parts at 
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the Keplerian distances of about 12 and 5 R©. The period 
variability may be interpreted as a spiralling-in or inward 
drifts in the innermost disc. The shortest observed period 
of 2.2 ±0.1 d may define the inner edge of the accretion 
disc at 4.8 ± 0.2 R ^, which agrees with the estimate of 
iMalbet et ail (|2005h from the interferometric observations. 
Within the disc-locking mechanism picture of the accretion 
processes, the rotational period of the star would be then 
close to 2 days. Although this was the interpretation of 
such oscillatory f eatures for TW Hya (|Rucinski et alj [20081 ; 
ISiwak et al-bonT ). Herbst (priv. comm.) argued that for the 
case of TW Hya, hot condensation in the disc could not eas- 
ily produce up to 20% flux variations observed in the visual 
spectrum which is dominated by the central star. This is in- 
deed a valid point, but we note that for the case of FU Ori 
the accretion disc luminosity overwhelms that of the central 
star by about 100 x. 

Two colour Stromgren v and b ground-based observa- 
tions of FU Ori confirmed the A(B — V) versus AV rela- 
tion obtained bv lKenvon et al l (|2000h . However, we obtained 
slightly redder colour index of the dominant variations, what 
may be due to a more outward hotspot location (causing 
the 8 d periodicity) during our observations. This would 
agree with the interpretation of the wavelength- amplitude re- 
lationship through different locations of the dominant vari- 
able flux, with longer periods produced by more external 
and cooler parts of the accretion disc. In generally this ef- 
fect may manifest itself within the wide-band MOST light 
curve and in the wavelet transform of the star variability 
as the amplitudes of the two drifting, oscillatory features 
markedly increased as their periods decreased. Future simul- 
taneous MOST and ground-based, multi-band observations 
conducted over a few months may help in investigation of 
the wavelength- amplitude relation for a range of possible pe- 
riods. Such observations should be supported by spectral 
synthesis disc models, to localize the dominant light sources 
at various inner disc concentric annuli. 

Z CMa, the complex binary consisting of a FUor and 
Herbig Be stars, was observed by MOST during 13 days 
close in time to the outburst maximum of the Herbig Be 
component. This led to the suppression of the FUor vari- 
ability component and dominance of the very chaotic spec- 
trum of oscillations of the Be component. The variability 
spectrum is somewhat similar to that of the Herbig Ae star 
HD 37806 (jRucinski et al.l I201Q ) . Future observations dur- 
ing quiescence level could help to disentangle the two light 
contributions. 
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